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glement in the composite state. Using a general quantum dynamical system we obtain 
\ the analytical solution when the atom initially starts from its excited state and the field 

in different initial states. Different features of the entanglement are investigated when the 
field is initially assumed to be in a coherent state, an even coherent state (Schrodinger 
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cate state) and a statistical mixture of coherent states. Our results show that the set- 
ting of the initial state and the Stark shift play important role in the evolution of the 
■ sub-entropies and entanglement. 

^. 

chaos and complexity and optical spatio-temporal dynamics}, 03.65.Ud {Entanglement 

+— > ' 

and quantum nonlocality (e.g. EPR paradox, Bell ?s inequalities, GHZ states, etc.)}, 
^ | 03.67.Hk {Quantum communication} 

> 

1 Introduction 

The quantum entropy and entanglement is a vital feature of quantum information. It 
has important applications for quantum communication and quantum computation, for 
example, quantum teleportation, massive parallelism of quantum computation and quan- 
tum cryptographic schemes [1-3]. Therefore, it is very essential and interesting how to 
measure the entanglement of quantum states. Recently much attention has been focused 
on the entanglement of the field and atom when the system starts from a pure state [4-16]. 
Also, in the context of the initial mixed state some studies have been reported [17-20]. In 
this context it was shown that calculating the partial entropies of the field or the atom can 
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be used as an operational measure of the entanglement degree of the generated quantum 
state. One finds that the higher the entropy, the greater the entanglement. Starting from 
an initial atom-field product state one can find perfectly entangled states between field 
and atom at certain later times even for initial coherent states with large photon number 
[4-8]. However, the time evolution of the field (atomic) entropy reflects the time evolution 
of the degree of entanglement only if one deals with a pure state of the system with zero 
total entropy. 

The quest for proper entanglement measures has received much attention in recent 
years [1,9]. From the identification and study of properties of such measures a gain of 
insight into the nature of entanglement is expected. In turn, their computation for par- 
ticular states provide us with an account of the resources present in those states. Because 
one needs to understand the best way to benchmark states for quantum information pro- 
tocols, here we examine the quantum partial entropies in the atom-field interaction for 
more general entangled states. From a practical point of view, an implementation of the 
quantum entropy will be used to measure the entanglement degree when the atom is as- 
sumed to be in its excited state and the field initially is in a coherent state, superposition 
state and a statistical mixture of two coherent states. As far as we are aware in the pre- 
vious investigations, that have dealt with the present problem, the initial system density 
matrix is taken to be a product of two states of the factored form. The atom is often 
taken to be in the excited pure state or mixed state and the radiation field is taken to 
being a pure state density matrix. So, the present task is a nontrivial issue, since we look 
at the mixed state entanglement from other direction taking into account the entropy of 
the field not equal to the entropy of the atom, in this case. To overcome such a difficulty, 
we employ the quantum von-Neumann entropy to measure the entropy of the atom while 
a numerical method will be used to calculate the quantum entropy of the field. 

The material of this paper is arranged as follows. In section 2, we find the exact 
solution of the system and write the expressions for the final state vector at any time 
t > 0. We investigate the quantum field (atomic) entropy and the atom-field entanglement 
in section 3. Finally, numerical results and conclusions are provided in section 4. 

2 The model 

The system we will consider here consists of a two-level atom interacting with a single- 
mode quantized field via k-quanta processes. The Hamiltonian in the rotating wave 



2 



approximation [8,20], can be written as (H — 1): 

H = H A +H F +H in , (1) 

where 

Hp = uja^a, 

H m = a^a((3 1 \g)(g\+(3 2 \e)(e\) + X(a^a_ + a k a + ), (2) 

where uj is the field frequency and u is the transition frequency between the excited 
and ground states of the atom. We denote by a and a* the annihilation and the cre- 
ation operators of the cavity field respectively. (3 1 and f3 2 are parameters describing the 
intensity-dependent Stark shifts of the two levels that are due to the virtual transition 
to the intermediate relay level, A is the effective coupling constant, a z is the population 
inversion operator, and a± are the "spin flip" operators, with the detuning parameter 
A = u — ku. 

Let us consider, the atom starting in its excited state |e), i. e., 

pg = |e><e|, (3) 

and we are going to assume that the initial single mode electromagnetic field inside the 
cavity is in a superposition state of the kind: 

Po = — + r 2 \ — a)(—a\ + r\a)(—a\ + r\ — a)(a\j , (4) 

where A = (1 + r 2 + 2r exp(— 2a 2 )) ,with a real. The parameter r takes the values —1, 
and 1, which corresponds to an odd coherent state, a coherent state and an even coherent 
state, respectively. While for certain classes of states a superpositions of coherent states, 
methods solely based on linear optical elements like beam splitters and photodetections 
could be found [21], an implementation covering other classes of entangled states remains 
a challenge. 

Also we want to see how different would the behavior of the system be if the input 
state is a statistical mixture of states \a) and | — a), i.e., 

f 1 

po = 2 (HH + 1 -«)(-«!)• (5) 

It is to be noted that when we put r = in equation (4) we get the same result as in Ref. 
13. It is expedient to expand the atom-field state in terms of the dressed states: 

|^/+' ) ) = sin# n |n, e) + cos# n |n + k, g), 

l^-'*) = cos# n |n, e) — sin 9 n \n + k, g), (6) 
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which are the eigenstates of the interaction Hamiltonian, where 



H\s,g) = E Q \s,g), 0<s<k 

H\^) = E^l^), (7) 



with the eigenvalues E and E±^ 



EF = u;\ [ n+-j+^ + +-[n(3 2 + (3 1 (n + k)]±Li n , 
E = (^i-y). (8) 



where 



n 



v n = — + -(/3 2 n- ^(n + k)), 

r. = (9) 



fi n is a modified Rabi frequency. The angle 6 n is given by 



e n = sin- 1 ( \ ] • (10) 



The unitary operator U t can be written as 

oo 

U t = £{exp(-^ ) )|^ n) )(^ ) | + exp(- ? tEi n) )|^ ( ! i) )(^ ( " ) |} 

n=0 
fc-1 

+ $>xp(-it£ )M>(#,s|. (11) 

s=0 

Despite being straightforwardly solvable in this way, the JC-model is well-known for the 
fact that the time-evolution of most expectation values is usually expressible only in series 
form. Having obtained the explicit forms of the unitary operator U t , for the system under 
consideration then the eigenvalues and the eigenfunctions can be used to discuss many 
features concerning the field or the atom. 

Bearing these facts in mind we find that the evolution operator U t takes the next from 



Pt 



Pi P2 
P3 Pi 
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where ( Pi ) nm = (n\ Pi \m), i = 1,2,3,4. ( Pl ) nm = A n (t)A* m (t), (p 2 ) nm = A n (t)B* rn _ k (t), 
(Ps)nm = B n -k(t)A*Jt), and (p 4 ) nm = B n _ k (t)B* m _ k (t). The coefficients A n (t) and B n (t) 
are given by 

a I \ ^ r -a r / \1 ( \ ■ sin XtfJ, \ 

A n {t) = q n C n exp[-iXtd + (n)\ \cos\tji n - in r 



Pn 



B n (t) = -iq n C n is n exp[-i\t5 + (n)] m ' '' 1 ' 



5+(n 



2^[n± R 2 (n + k)], when R ^ 



where 



when (3 i = 0, 



Cn = [-^(l + r(-ir)], 



(12) 



for the initial condition (3), while for the initial condition (4) is 

C n = [±={5 l + {-\Y5 )l (13) 

with <5j, 5j satisfying the two following condition, (a) 5i = 5j = (5,j) 2 = (5j) 2 = 1, and (b) 
Si.Sj = 0, 

00 00 n 

l«) = E^|n) = Ee" Q2/2 ^|n). (14) 

n=0 n=0 Vnl 

With the final state obtained, any property related to the atom or the field can be calcu- 
lated. Employing the reduced density operator for the atom or the field, we investigate 
the properties of the entropies (S a , Sf) and hence entanglement. 



3 Entropy and subentropy 

There is growing interest in the roles of nonadditive measures in quantum information 
theory. Inadequacy of the additive Shannon von-Neumann entropy as a measure of the 
information content of a quantum system has been pointed out [22]. Also, there is a the- 
oretical observation [23] that the measure of quantum entanglement may not be additive. 
Despite the fact that the basic idea of quantum entanglement was acknowledged almost 
as soon as quantum theory was discovered, it is only in the last few years, that considera- 
tion has been given to finding mathematical methods to generally quantify entanglement. 
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In the case of a pure quantum state of two subsystems, a number of widely accepted 
measures of entanglement are known. However, the question of quantifying the degree of 
entanglement for general mixed states is still under discussion. Let us now briefly repeat 
some of the key underlying definitions. The entropy S of a quantum-mechanical system 
described by the density operator p is defined as follows: 



S = -Tr{p lnp}, 



(15) 



where we have set the Boltzmann constant K equal to unity. If p describes a pure state, 
then S = 0, and if p describes a mixed state, then S 7^ 0. Entropies of the atomic and 
field sub-systems are defined by the corresponding reduced density operators: 

Sa( f ) = -Tr aU) {p aU) lnp a{f) }. (16) 

Taking the partial trace over the field, the reduced atomic matrix can be written as 



Pt 



Tr f {p) 



Pee Peg 
Pge Pgg 



E \Mt)f 

n=0 



E A n (t)B* n+k (t) 

n=0 



E B n+k (t)A* n (t) E \B n+k (t)\ 2 

\ n=0 n=0 



(17) 



Thus we rigorously obtain the quantum atomic entropy in the following form 

S{p a t ) = -A«(t)logA«(t) - X a _(t)\og\ a _(t), (18) 

where A"(t) is given by 

A° ± (0 = I {l ± ^(2Pee(t)-l) 2 +MPe 9 (t)\ 2 } " ( 19 ) 

In this case, the probability of finding the atom in its excited or ground states are expressed 
as the diagonal element of the reduced atomic density matrix, i.e., 

Pu(t) = (i\Pt\i), i = e,9 (20) 



and the off-diagonal element p eg (t) is given by 

PuM = 



i = e,g. 



(21) 



Taking the partial trace over the atomic system, we obtain the reduced density operator 
in the form 

p{ = tr A p(t), (22) 
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with its (pt)nm element given by 

(p{)nm = A n (t)A* m (t) + B n _ k (t)B* m _ k (t). (23) 

From this equation, it is difficult to obtain the eigenvalues of the reduced density operator 
for the field, in this paper we will evaluate them numerically. 

4 Results and conclusion 

We study the temporal behavior of the atom-field system in the JC-model for the cavity- 
field prepared initially in different forms. As an example we may consider a simple initial 
condition for the atom to be in the excited state and the field in a coherent state or 
a superposition of the coherent state or in a statistical mixture of two coherent states 
(equations (3) and (4)). Among the family of mixed quantum mechanical states, special 
status should be accorded to those for a given value of the entropy and have the largest 
possible degree of entanglement. The reason for this is that such states can be regarded 
as mixed-state generalizations of Bell states, the latter being known to be the maximally 
entangled two-qubit pure states. Hence, this kind of mixed states could be expected to 
provide useful resources for quantum information processing. At this end, we have the 
plot of the quantum partial entropies (S a , Sf) relative to these different initial states of 
the atom-field, as a function of the scaled time Xt/n, taking into account the two-photon 
process (k — 2), in order to investigate the Stark shift effects. 

We assume a fixed value of the initial mean number of quanta n — 16 and different 
values of Stark shift parameter R (namely, /3 1 — f3 2 — 0, i.e. in the absence of Stark shift 
in figure 1, R = 0.5 in figure 2 and R = 0.3 for figure 3). Furthermore, the detuning 
parameter is taken to be zero, and in figure la we set r = (coherent state), figure 16, we 
set r = 1 (even coherent state) and figure lc (a statistical mixture state). In the absence 
of the Stark shift, it is observed that the quantum field entropy and the quantum atomic 
entropy have the same values due to the initial coherent state r = 0, (see figure la). 
This behavior is similar to that obtained in the standard two-photon two-level systems 
obtained previously (see for example [4-5]). It is observed that the entropy evolves with a 
period tt/X, when t = mr/A, n — 0, 1, 2, ...), the quantum field entropy evolves to zero and 
the field is completely disentangled from the atom, while for t — (n+ |)vr/A, it evolves to 
the maximum value, and the field is strongly entangled with the atom. 

The situation is completely changed when we consider an even coherent state i.e r — 1. 
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Figure 1: The entropy for the atom S a (t) (solid line) and the field entropy Sf(t) (dashed 
line) as a function of the scaled time Xt/ir of the particle initially prepared the excited 
state and the field initially prepared in: (a) a superposition state SS (r = (coherent 
state)), (6) a superposition state (r = 1 (even coherent state)) and (c) a statistical mixture 
(MS) of coherent states | a) and | —a) (n = 16). 
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Although the quantum field entropy is still equal the atomic entropy, the quantum entropy 
in this case has minimum values at Xt = ™, (n = 0, 1, 2, ...) i.e at half of the revival time, 
also instead of the two pre-minimum values observed in figure la we have here only 
one pre-minimum value between each two consecutive minima. This effect is due to the 
interference between the two coherent states in the superposition, and can be understood 
by looking at the photon number distribution of the initial field. Hence this signal gives 
a clear measure of the remaining degree of coherence between the two components of the 
Schrodinger cat state, while the signals present in both cases are due to intrinsic revivals 
of each component individually. Because of this enhancement it is possible to have the 
generation of well-defined Schrodinger cat-like states during the evolution of the field in 
the two-photon process case model [24]. We would like to remark that the approach to a 
pure state at half of the revival time occurs in the ordinary JC-model [25], if we start with 
the field in a pure state. Let us now come to specific numerical examples to investigate the 
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Figure 2: The same as in figure 1 but with the Stark shift parameter R = 0.5 

influence of the statistical mixture on the evolution of the quantum field entropy and the 
quantum atomic entropy (see figure lc). An understanding of interaction between an atom 
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and an electromagnetic field has been possible in recent years through the introduction of 
the statistical mixture state picture [26]. By looking to the statistical mixture state one 
has a clear physical understanding of what are the parameters involved in such expression 
and what is going to neglect in order to go from a pure state to a mixed state. On the 
other words, the state of the initial field is an equally-weighted statistical mixture of two 
coherent states, which is a special class of the Schrodinger cat state. In this it has 

been already discussed [20], the quantum field entropy Sf(t) is greater than the quantum 
atomic entropy S a (t) (see figure lc), S a (t) reach its maximum values at half of the revival 
time, while Sf(t) evolves to minimum values. We note a little deviation from ordinary 
Rabi oscillations, due to the statistical mixture case. It is seen that entanglement evolves 
depending on the initial preparation of atoms, however from figures la and lc we see that 
the quantum atomic entropy has similar behavior in both superposition and statistical 
mixture of coherent state. 
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Figure 3: The same as in figure 1 but with the Stark shift parameter R = 0.3 
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For comparison purposes, we have chosen to set some different values of the Stark 
shift parameter R and the other parameters are the same as in figure 1. The outcome is 
presented in figure 2 (where R = 0.5) and figure 3 (where R = 0.3). We note a stronger 
modulation in the oscillations and a clear departure from ordinary Rabi oscillations being 
verified as the Stark shift parameter takes values far from the unity (see figures 2 and 
3). Further studies about such a situation have been carried out and are not displayed 
here. It is interesting to refer here to the fact that the Stark shift creates an effective 
intensity dependent detuning An = (3 2 — Pi [27]. When f3 2 — Pi, Ajv = 0, in this 
case, the Stark shift does not affect the time evolution of the quantum entropy. As is 
visible from the figures, the effects of the dynamic Stark shift are more pronounced when 
R deviates from unity. Interestingly, when R is decreases, the values of the maximum 
entropy are decreased. Periodic models therefore may be more robust in this sense. Also, 
with decreasing the parameter R, the evolution period of the entropy as well as the 
subentropies is decreases (see figure 3). The sensitivity becomes even more clearly visible 
when we take small velues of the Strak shift parameter. It is worth mentioning that the 
Strak shift effect has the same impact for both the field entropy and the atomic entropy. 

When the atomic entropy is calculated, we note that the terms involved are of the forms 
(U(t)a\U(t)a) , {U(t)a\ — U(t)a) and | (U(t)a\U(t)a) | 2 for the case of the mixture. These 
terms in particular do not differs from those of the case of the coherent state which has 
been discussed earlier. Therefore the temporal evolution of the entropy for the atomic 
system alone in the case of the Schrodinger cat sate (S) mimics the evolution of the 
entropy in the pure coherent state as can be seen from comparing figures (lc) and (la). 
However when we consider the entropy for the field we note that terms of the form 
{±U(t)a\ ± U(t)a) with all combinations. These terms are the ones that appear in the 
case of the superposition of the two states (r = 1 in equation 4). Hence the resemblance 
between the figures for the entropy of the field in the mixed state of figure (lc) and the 
case of the initial superposed states of figure (lb). This may demonstrate relevance of 
investigating the entropies of the subsystems and their relation to entanglement. 

In summary, we have shown in this paper that the final analytical expression of the 
composite density matrix along with its overlap matrix elements can be used to obtain 
the quantum field entropy and the quantum atomic entropy. This is accomplished by 
choosing to study the system in the representation in which the marginal initial density 
matrices are assumed to be in a coherent state, superposition states and statistical mixture 
states of two coherent states. We present different numerical examples to elucidate the 
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effects of these different settings. Explicit computations are presented for different values 
of the Stark shift parameter. Our results show that the superposition of coherent states 
and Stark shift play an important role in the evolution of the quantum entropies in the 
two-quanta JC-model. In the coherent state, the quantum field (atom) entropy reaches 
its maximum values at the half of the revival time and the field and the atom are strongly 
entangled, while in the even coherent state, the entropies at the same time evolve to 
the minimum values (zero) and the field and the atom are strongly disentangled. In the 
statistical mixture state, (S a 7^ Sf), the entropy for the atom reach to the maximum 
values at half of the revival time, while Sf evolves to its minimum values. The significant 
effect of the Stark shift parameter appears when R deviates from unity. The more R 
deviates the more the two systems are weakly entangled. 
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